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We show that the 511 keV gamma ray excess observed by INTEGRAL/SPI can be more robustly 
explained by exciting dark matter (DM) at the center of the galaxy, if there is a peculiar spectrum 
of DM states xo, Xi an d X2, with masses Mo ~ 500 GeV, Mi <, Mo + 2m e , and M2 = Mi + 
5M £ Mo + 2m e . The small mass splitting 5M should be £ 100 keV. In addition, we require at 
least two new gauge bosons (preferably three), with masses ~100 MeV. With this spectrum, \i is 
stable, but can be excited to \2 by low-velocity DM scatterings near the galactic center, which are 
Sommerfeld-enhanced by two of the 100 MeV gauge boson exchanges. The excited state X2 decays 
to xo and nonrelativistic e + e~, mediated by the third gauge boson, which mixes with the photon 
and Z. Although such a small 100 keV splitting has been independently proposed for explaining 
the DAMA annual modulation through the inelastic DM mechanism, the need for stability of Xi 
(and hence seqestering it from the Standard Model) implies that our scenario cannot account for 
the DAMA signal. It can however address the PAMELA/ATIC positron excess via DM annihilation 
in the galaxy, and it offers the possibility of a sharper feature in the ATIC spectrum relative to 
previously proposed models. The data are consistent with three new gauge bosons, whose couplings 
fit naturally into a broken SU(2) gauge theory where the DM is a triplet of the SU(2). We propose 
a simple model in which the SU(2) is broken by new Higgs triplet and 5-plet VEV's, giving rise 
to the right spectrum of DM, and mixing of one of the new gauge bosons with the photon and Z 
boson. A coupling of the DM to a heavy Z' may also be necessary to get the right relic density and 
PAMELA/ATIC signals. 

PACS numbers: 98.80.Cq, 98.70.Rc, 95.35.+d, 12.60Cn 



I. INTRODUCTION 

Dramatic developments in observational astronomy 
have started to alter our picture of dark matter (DM); in- 
stead of being a single state, observations have indirectly 
suggested that DM could be a multiplet with small mass 
splittings. In ref. [l| it was argued that such a scenario 
can explain not only the positron/electron excess recently 
indicated in the 10-100 GeV region by the PAMELA 
and (to some extent) HEAT 0] experiments, and a simi- 
lar one in the 500-800 GeV region seen by ATIC 0] and 
PPB-BETS d, but also the 511 keV gamma rays ob- 
served by INTEGRAL/SPI @, Q, the WMAP haze 
and the annual modulation observed by DAM A/LIBRA 
Q. The connections were further explored in [lfj. The 
common link between them is that they can be explained 
in terms of DM undergoing enhanced scattering and sub- 
sequent annihilation into light bosons which decay to 
e + e~ 111 ], or else exciting a DM state with a small mass 



splitting above the ground state, which might decay back 
to the ground state and e + e~ if the splitting is greater 
than 2m e . In the case of DAMA, the mass splitting pro- 
vides the kinematics which would enable DAMA to be 
more sensitive than other experiments to a necessarily in- 
elastic collision. The ideas of excited dark matter (XDM) 
[l2| and inelastic dark matter (IDM) [l3| were proposed 
before ref. but the latter took the step of trying to 
unify them into an appealing theoretical framework, and 
to use it to also explain the excess positron/electron re- 
sults of PAMELA/ATIC. 



The unified description of dark matter has one short- 
coming, however. In ref. [3], it was shown that the XDM 
mechanism falls short of being able to reproduce the ex- 
perimental observation by nearly three orders of magni- 
tude, even if the galactic DM scattering cross section <j ga \ 
saturates the unitarity limit in the s-wave contribution. 
The main loophole for circumventing this conclusion was 
to hope that higher-Z partial waves could increase <7 ga i 
by a factor of 300. Such an enhancement was argued to 
be unlikely in ref. [l4j]> and we will show that this argu- 
ment is borne out in the class of models proposed by ref. 
[l[. Thus the XDM explanation of the 511 keV excess 
remains unrealized. 1 

We note that the 511 keV anomaly is not just the find- 
ing of the INTEGRAL experiment, but it was first ob- 
served in 1972, and has been seen in four subsequent 
balloon- and satellite-borne experiments The obser- 
vation is thus quite credible, and so far is lacking any 



1 This conclusion depends on the small error bars for the necessary 
value of Ugai estimated by ref. [l7ll . However there seems to be a 
larger uncertainty in the DM density at the center of the galaxy, 
n c . Since positron production in the XDM scenario scales like 
v? c , the original XDM idea with SM = m e might be salvaged 
if n c is greater by a factor of 17 (even more if the unitarity 
bound is not saturated) than in the model found to be preferred 
in the best fit to the INTEGRAL data by ref. [JJ. It would 
be worthwhile, though beyond the scope of the present work, to 
further investigate this point. 
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FIG. 1: Left: DM spectrum needed in present work to ac- 
count for the INTEGRAL/SPI observations; right: spectrum 
suggested by ref. 



highly convincing astrophysical or particle physics expla- 
nation, although attempts have been made using positron 
emission from low-mass x-ray binaries (LMXB's) [15j, or 
annihilation of light MeV-scale dark matter (see for ex- 
ample [H, [T3|) as well as decaying relics 0, [U| . The 
LMXB hypothesis is predicated on a supposed correla- 
tion between the asymmetry in the disk component of 
the 511 keV gamma rays and the distribution of bright 
LMXB's, but this has been criticized on several grounds 
in ref. [19| (although the discrepancies might be amelio- 
rated if positrons produced in the disk can be transported 
to the bulge before annihilating [20j]). Therefore it is still 
interesting to find a technically natural particle physics 
explanation for the 511 keV line, regardless of the other 
experimental anomalies. This was the primary motiva- 
tion for the present work. The fact that our positive find- 
ing for the INTEGRAL anomaly is consistent with the 
general framework outlined in ref. [l[ , for also explaining 
the other experiments (apart from DAMA), heightens its 
interest. 

Our new twist for making XDM viable is to have a mass 
splitting SM which is much smaller than m e between the 
middle DM state and the heaviest one, as illustrated in 
figure [TJ These are assumed to be Majorana fermions, 
Xi- We will design the model so that xi cannot decay 
directly to the ground state xo! thus both xo and Xi 
are stable, xi undergoes Sommerfeld-enhanced scatter- 
ing P, [U [22|, [H[ at the galactic center (where the DM 
density is highest) through multiple exchange of a light 
gauge boson with off-diagonal coupling X^T^-B^Xi, 
as envisioned in ref. [l|. The produced xi states subse- 
quently decay through another vector B'^ which mixes 
with the photon to produce e + e~. These processes are 
shown in figure O The key ingredient which makes the 
scattering efficient enough is the small mass splitting SM. 
Previous attempts to implement the XDM mechanism 
have failed because a larger splitting ~ m e was assumed, 
and this makes the excitation rate too small. The figure 
of merit is the ratio of the predicted rate of e + produc- 
tion to that observed, in the Ith partial wave, assuming 
the unitarity bound is saturated [14j : 
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FIG. 2: Left: Sommerfeld-enhanced scattering X1X1 — * X2X2- 
Right: decay \2 — ► Xoe + e_. 



It is suppressed by the Maxwell-Boltzmann distribution 
at the threshold for production of the excited state. Eq. 
(fTj) was predicated on the assumption that the kinetic 
energy of each DM particle must be sufficient to produce 
one electron or positron. Since the characteristic DM ve- 
locity vq is fixed, M must be sufficiently large to avoid 
the Boltzmann suppression, but the rate also scales like 
M ~ 4 . Even with the optimal value of M = m e /2vQ, one 
would need to more than double the estimated value of 
vq = 6 x 10 _4 c to make Rq = 1. Such a large change 
seems to be well outside the range of uncertainty in 
the current understanding of the DM velocity distribu- 
tion. However in our scenario, m e is replaced by the 
smaller SM, which can significantly ease this tension. 
We will show that ^ Ri can be 1 as required without 
changing vq, and keeping M ~ 500 GeV as desired for 
PAMELA/ATIC, if SM < 100 keV. It is intriguing that 
the same splitting has been advocated previously to ac- 
count for the DAMA signal. However we will show that 
the need for stability of xi means that it cannot inter- 
act with baryons or leptons at detectable levels; thus our 
proposal does not seem to be compatible with the IDM 
explanation for DAMA. 

In the remainder of the paper we will give details of the 
computation of the excitation rate (section II), tighten 
the case against the large mass gap SM = m e scenario in 
section III, and show how a smaller value SM ~ 100 
keV can improve the situation in section IV. In sec- 
tion V we build a simple particle physics model of DM 
which can accommodate our findings, and also address 
the PAMELA / ATIC /PPB-BETS observations. Section 
VI gives a brief account of the cosmological implications 
of the model. Conclusions are given in section VII. 



II. METHODOLOGY 

To obtain the desired result for the 511 keV signal, it 
is important to nearly saturate the unitarity bound in at 
least one partial wave (or to have significant scattering 
up to high values of I, but we will show that this does 
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not seem to be possible in the present context). The 
Sommcrfcld enhancement which can occur at low DM 
velocities is crucial for getting such strong scattering. We 
follow the quantum mechanical treatment of appendix 
A. 4 of ref. [1] to compute this effect. 

Because the gauge coupling gxiR^Xi is assumed to 
be off-diagonal, we have two states |1) = |xi>Xi); |2) = 
1X25X2)5 whose interaction Hamiltonian has the matrix 
form 



Vij = 







-ae 



-» r /r 2SM 



(2) 



where a = g 2 /4ir and u is the mass of B^. The wave func- 
tion for the two-state system (with components labeled 
by index i) in the CM frame is = J2i p i{ C0S 9) R ki( r )' 
where k is the initial momentum. Defining $/.i(r) = 
R l k i/r, the Schrodinger equation is 2 



1 a>" 4- ( l SL±As 



Vi. 



(3) 



The equation is solved by the shooting method, where 
<j> ; ^ r i+1 (*) near r = for some complex number b, 
which is then adjusted so that there are only outgoing 
and not incoming waves in <E>/2 as r — > 00. 

To extract the scattering amplitudes, we decompose 
the numerical solution into incoming and outgoing waves, 
<&n.> Q\ n 2 and $"2'' Partial wave unitarity implies the 
conservation of flux, fc^"! 2 = k\^f\ 2 + k'\^f\ 2 where 
k' 2 = k 2 — 2M%6M, which we use as a check on our 
numerics. The fraction of incoming |xi, Xi) states which 
gets converted to the |x2, X2) final state is thus 
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in the lib. partial wave. This must be integrated with the 
Maxwell-Boltzmann distribution Nv 2 e~ v / v ° (or some 
more sophisticated distribution function, as we discuss 
below) to find the thermally averaged cross section 
(cgai^rei)- Doing so modifies the unitarity-saturating es- 
timate ([1} to read 



Ri -> 4.5 x 10" 



21 + 1 /500 GeV 



Mi J J ut 



due u fi(u) 



(5) 

where u — v 2 /vq, vt — d2bMjM\ is the threshold ve- 
locity for X2 production and v csc (r) = 700 km/s is the 



escape velocity at r = 0.4 kpc, the outer edge of the re- 
gion where INTEGRAL sees excess 7 emission [l| ■ Since 
Vq = 180 km/s, the error in extending the upper limit of 
integration to 00 is small. 



III. CONSTRAINTS ON LARGE 5M MODEL 

We have manually scanned the parameter space of the 
model to try to maximize the fraction of excited state 
particles, //, first starting with the original class of mod- 



els with the "large" mass gap, SM 



The goal 



here was to see if any parameters could be found such 
that high partial waves could contribute, thus overcom- 
ing the small prefactor Ri = 3.4 x 10~ 3 (2Z + 1) at the 
optimal mass M — m e /2w§ and vq = 180 km/s We 
considered the limit in which the gauge boson mass 
H can be neglected. Otherwise the range of the interac- 
tion is reduced; this can only decrease the contribution 
from higher-? partial waves, which correspond to scatter- 
ing at large impact parameter. In that case, the relevant 
dimensionless parameter turns out to be 



M x a 2 
26M 



(6) 



(see footnote [2]) For T ^> 1, the s-wave dominates the 
cross-section, while for V <C 1, a range of partial waves 
contribute significantly. However, in the latter case even 
the largest contributions fall short of the needed unitarity 
limit, so the total cross section is not actually enhanced. 
Working at the optimal mass M — m e /2v 2 , we have 
computed 



I' = e Ut 



I 



(21 + 1) 



due u fi(u) 



(7) 



[see eq. ©J as a function of T. This quantity would take 
the limiting value + 1) f° r an the partial waves 

which reach the unitarity limit // = 1. At the optimal 
mass M = m e /2«Q, u t = 4, and ref. [14| shows that /' 
should have the value (3.4 x 10 -3 ) -1 = 300 in order to 
match the INTEGRAL observations. Fig. 3 shows that 
in the actual model, /' reaches a maximum value of 0.3 
near T = 0.5, far below what is needed. 

We have thus established that the high-? loophole for 
the large 5M scenario does not work. Another way of 
enhancing the effect would be to take advantage of its 
strong dependence on ^o- One would need to boost vq by 
the factor (0.3 x 3.4 x 10" 3 )" 1 / 7 ^ 2.7, giving v = 480 
km/s instead of 180 km/s. On the other hand, more 



2 For the numerical solution it is useful to rescale r = (a/28M)x 
and define the dimensionless variables T = Mi a 2 / '(2SM), T = 
(ak/2SM) 2 , 77 = a[i/2SM, so that the Schrodinger equation 
takes the form + [1(1 + l)/x 2 + TV]* = T*, with the di- 



mensionless potential V = ( 







). We must have 



T > r for the initial state to have enough energy to produce the 
heavier |X2,X2) final state. 



3 In terms of ref. [J , this actually requires looking at the scattering 
X0X1 — > X1X2 since X1X1 — * X2X2 would have SM 2in e and 
be even more suppressed. Ref. 11411 has SM m e by virtue of 
charged intermediate states, which allows each incoming \ to be 
excited by only m e rather than 2m e . 



FIG. 3: Left: Enhancement factor I', eq. 0, as a func- 
tion of r = Mia 2 /25M , in the disfavored case SM = m e , 
where the new gauge boson is taken to be massless and M 
is optimized. The maximum enhancement is too small to 
match observations. Right: integrands for successive partial 
waves contributing to /', versus u/a, where u = v 2 /vq and 
a — 26M/Mvq. One would need significant contributions 
from many more partial waves to make I' large enough. 



sophisticated estimates of the dark matter distribution 
indicate that v is smaller than the fiducial value, rather 
than larger, in the center of the galaxy [TJ, [r7| . Dif- 
ferent models of the DM distribution function n(r, v) 
give radially-dependent average velocities vo(r) which de- 
crease toward r = 0. Therefore using the constant value 
for wo which best describes the bulk of the galaxy al- 
ready overestimates the efficiency of DM excitation near 
the center, and the large SM possibility seems to be ruled 
out (see however footnote [T]). 



IV. THE CASE OF SMALL SM 

Now we turn to the main point, that smaller values of 
the DM mass splitting SM ~ 100 keV can overcome the 
problem of too small a signal, without any need for in- 
creasing the DM velocity vq. Redoing the analysis of ref. 
[3] for general SM, one finds that for the optimal DM 
mass, Ri is enhanced by an extra factor of (m e / SM) A . 
We can therefore achieve the desired effect if unitarity 
is nearly saturated only in the s-wave (or other low-/ 
contributions), with SM = 0.24 m e = 120 keV. This es- 
timate applies for the optimal mass M = SM/2vq = 170 
GeV. However we can make the mechanism work at larger 
M, as desired for getting the unified explanation of the 
PAMELA / ATIC /PPB-BETS observations Q, by mak- 
ing SM only moderately smaller, as we will now show. 

We have done a preliminary exploration of the param- 
eter space, to see what can be achieved in the concrete 
framework at hand. We defer a more comprehensive anal- 
ysis to the future; here we will just present a working 
example. To get a large enough effect, it is important to 
vary the mass fj, of the exchanged gauge boson B^. Phys- 
ically, this is due to resonant scatterin g w hen a bound 
state of nearly zero energy forms [l], l2lT |23| . This effect 
is only possible for a finite-range potential such as the 
Yukawa type. Generally, we find enhanced scattering for 
larger values of Y = Mia 2 /26M, which is not surprising 
since the interaction strength is a, and there is an optimal 
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FIG. 4: Enhancement factor / = e~ ut I' — e~ a I' versus a = 
2SM/(Mv 2 ) for the small SM model, with T = 10, T] = 1.2. 
Successive contributions from partial waves Z = — 10 are 
shown. I = 1.3 (dotted line) is the value indicated by the 
INTEGRAL observations. 

(though not sharply peaked) value of the dimensionless 
parameter 

r) = afj,/25M (8) 

giving the resonant effect. 

As an example, we present the case where T = 10, 
r] = 1.2. We consider the quantity I = e~ Ut I' = e~ a I' 
rather than V of eq. (7}, where a = 25M/(Mv%). In the 
large SM scenario, the best one could do was to optimize 
M such that a = 4, but for small SM, we can obtain much 
larger results at smaller values of a. In contrast to the 
large SM case, we can hold M fixed at a a value which is 
larger than optimum explaining for the 511 keV line, but 
more interesting for simultaneously explaining the other 
potentially DM-related anomalies. From eq. it follows 
that at M = 500 GeV and vq/c = 0.0006, one only needs 
I = 1.3 to explain the INTEGRAL observations. 

The enhancement factor / is plotted as a function a in 
figure 21 where the convergence of the successive partial 
wave contributions is shown. The needed value of I = 
1.3 can be obtained for a = 0.95. For M = 500 GeV 
and vq/c = 0.0006, this implies the mass splitting SM = 
86 keV. From T = 10 and 77 = 1.2, we infer that the 
fine-structure constant of the new gauge coupling is a = 
0.0017, and the mass of the gauge boson is fi = 120 
MeV. 

So far, we have assumed that the initial state xi nas a 
particular number density, namely that which was used 
in the analysis of 17] to determine the cross section 
CTgai = 10" 28 (Mi/TeV) 2 cm 2 needed to explain the IN- 
TEGRAL anomaly. However, with three nearly degener- 
ate DM states, one expects the xi density to be 1/3 this 
value, given that X2 decayed into xo an d not xi- The sig- 
nal is proportional to the integral of (cr ga iv re i)?T. 2 (r) along 
our line of sight, so this would require the cross section 
to be nine times higher than we have assumed. The ad- 
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ditional suppression can be counteracted if the DM den- 
sity in the central region r < 0.4 kpc is three times higher 
than assumed in the analysis of ref . [llj . Even apart from 
any uncertainties in the shape of n(r), we note that it is 
normalized to the local energy density n(8.5 kpc) = 0.3 
GeV/cm 3 , which is estimated to be uncertain by a fac- 
tor of 2 in the upward direction [24j . If this uncertainty 
works in our favor, we only need an additional factor of 
1.5 enhancement in the central region. This is a modest 
shortfall, since in the four models considered by ref. [TtJ , 
n(0.4 kpc) varies by a factor of 27. Thus the smaller xi 
density does not seem to pose a serious problem for the 
model. On the other hand, the missing factor of 300 in 
the large SM case looks more daunting. 

We have also neglected the probable r-dcpcndcnce of 
vo(r) discussed at the end of the previous section, which 
also tends to reduce the predicted signal. There is not yet 
a consensus on the precise form of vo(r), but if it proves 
to give a significant reduction, this can be compensated 
to some extent by taking more optimal (smaller) values 
of the DM mass. It is also possible that more favorable 
examples at large M\ exist, which would make it inter- 
esting to perform a wider and more systematic search of 
the model parameter space (r, rj) than we have been able 
to do so far. 
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FIG. 5: Energy spectrum of e 4 
circles), taken from ref. [3]. 
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data. We will come back to this issue in the next section. 

Let us show how a simple SU(2) model could account 
for the observations. If Xa is a real triplet, its gauge 
interactions are 



ge^XaB^^Xc 



(9) 



V. ATIC AND MODEL BUILDING 

The simplest and most theoretically appealing way of 
getting three DM states Xi an d several gauge bosons is 
to assume that the gauge symmetry is SU(2) and that Xi 
transforms as a triplet. We will now explore the conse- 
quences of this hypothesis, with particular attention to 
its implications for the ATIC/PPB-BETS excess e+e~ at 
300-800 GeV. Of course we will also maintain all ingredi- 
ents needed for our successful explanation of the 511 keV 
anomaly. 

Already to explain INTEGRAL, we required two gauge 
bosons, B which couples to X1X2 to facilitate XiXi ~~ > 
X2X2, and B' which mixes with the photon and medi- 
ates the decay xi — * Xa e+e ~ ■ However, with only these 
states, it is difficult to produce a rather sharply-peaked 
excess of high-energy leptons on top of a lower energy 
continuum, which is suggested by the ATIC data and 
which we have reproduced in figure O (The PPB-BETS 
data are consistent, but not so strongly suggestive of the 
peak.) The process specified by ref. [l[ in this regard is 
shown in fig. [51 where the final state B' bosons subse- 
quently decay to e + e~. The annihilation initially pro- 
duces back-to-back B's, each carrying energy Mi, which 
should give a leptonic spectrum that is rather uniformly 
distributed in energy. It would be desirable to have an ad- 
ditional channel which produces a single nearly monoen- 
ergetic e + e _ pair, each lepton having energy Mi. Such a 
spike would be significantly broadened by Coulomb scat- 
tering of the primary particles in the galactic medium 
J25|-[28j], possibly giving the peak-like shape in the ATIC 



Letting a, b, c = 0,1,2, then in the previous notation, 
B = B, Bi = B', B 2 = B". A simple way to get 
kinetic mixing of the photon with B' is by including a 
Higgs triplet A a in the dark sector, and the dimension-5 
operator 



_Y> iV B a A 



(10) 



where is the field strength of the SM weak hyper- 
charge gauge field. If we assume that only Ai gets a 
VEV, this generates the required mixing of B' with the 
photon and the Z boson. It is straightforward to work 
out the transformation that diagonalizes the kinetic term 




FIG. 6: Sommerfeld-enhanced annihilation channel for ex- 
plaining the PAMELA/ ATIC/PPB-BETS anomalies. Note 
that B' will decay into e + e~ . The process suggested by ref. 
[1( is shown. 



of the B', A (photon) and Z boson. If // is the mass of 
the B 1 , we find that the flavor states B',A,Z are related 
to the mass eigenstates B', Y, Z by 

A = A- ecos6 w B' + 0(e 2 ) 
B' = & + esm6 w Z + 0(e 2 ) 

u' 2 - 

Z = Z - esm9 w ^B' + Ote 2 ) (11) 
mj 

where e = 2Ai/A and 9w is the Weinberg angle. It is 
important for this class of models that B' does not mix 
with the photon A at O(e) [in fact neither does it mix 
at higher orders], because otherwise X2X0 would acquire 
a coupling to A which would be just as strong as the 
coupling of B' to e+e~. In that case, INTEGRAL would 
observe a narrow line from xi ~ ¥ Xo7 at energy SM ~ 1 
MeV, in addition to that at 511 keV, but this of course 
has not been seen. The A-B' mixing implies a decay 
rate of t _1 ~ a cm aMe 2 (m e / /x) for \i ~ * Xo e+ e~, and 
a corresponding lifetime of order 10~ 5 s if e ~ 10~ 4 . For 
e ~ 10~ 4 and (Ai) ~ 10 GeV as we will find below, the 
scale A is - 50 TcV. 

With the breaking pattern (A„) = <5 a iA needed in eq. 
(TIT)]) , the kinetic term (D^A) 2 of the triplet only gives 
mass terms for B and B" , 

g 2 [A a A a B b B b - (A a B a ) 2 } = g 2 A 2 (B 2 + B" 2 ) (12) 

(where g is the dark SU(2) gauge coupling), so we must 
also include symmetry breaking from another new Higgs 
field. One simple possibility is a symmetric traceless ten- 
sor, E ab , the 5D representation, whose VEV is only in 
the E02 = S20 or £00 = — ^22 components. In fact, a 
global rotation around the 1 direction (which is the sub- 
group of SU(2) left unbroken by our choice of (Ai) as 
the triplet VEV) can conveniently put the 5-plet VEV in 
the diagonal Son = —£22 components alone. The group 
generators in the 5D representation can be written as 

T ab,ce = iitadcSbe + ZbdeSac)- The VEV E = ~E 22 = S 

generates a mass term proportional to 

-g 2 B d B f (X ab )TX ce Ti <hi (X hi ) = g 2 E 2 (B 2 +B" 2 +2B' 2 ) 

(13) 

The resulting gauge boson mass spectrum is then /i = 
fx" = gVW+A 2 and p! = gV^E. To get the desired 
radiative mass splittings of the %'s below, we will need 
to assume that A < E, so that \i < yl '. 

An important property of the interaction © (and in- 
deed all the interactions in our model) is that it preserves 
the Z 2 symmetry B -> —B, B" -> —B", xi -> ~Xi 
which is needed to keep Xl stable, if we assume that 
Ao and A2 are also charged under the Z 2 - The dis- 
crete symmetry is thus unbroken by the VEV of Ai. 
Note that the nonabelian cubic and quartic interactions 
of the gauge bosons have the schematic form BB'B" and 
(B 2 + B" 2 )B' 2 , which also respect the Z 2 symmetry. A 
notable consequence of the symmetry is that B, B" do 
not acquire any couplings to SM matter which would 
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FIG. 7: Radiative corrections to the \ a masses, which leave 
Xo and X2 degenerate, fi and // are the masses of the B/B" 
and B' gauge bosons, respectively. 

allow them to decay into e + e~, nor to cause nuclear re- 
coil in direct DM detection experiments. Therefore even 
though the excitation xi X2 has the right kinematics 
for the IDM explanation of DAMA, the putatively ex- 
changed B-B" boson cannot interact with the detector 
(nor can xo ~ > Xi work, both for this reason and because 
of the larger mass splitting). The remaining excitation 
Xo — » X2 is possible from the point of view of the interac- 
tions of the exchanged B' boson, but the mass splitting 
~ 2m e is too large for it to proceed at a detectable rate. 

Next we consider the spectrum of the DM triplet. The 
tree level mass Mx a Xa gets split by ha((J,'—fJ,) [29| by the 
diagrams of fig.0 4 We are assuming that fi < p' since, as 
we will show, this is what gives the desired DM spectrum. 
Figure [7] indicates that the radiative correction gives xi a 
mass which is larger than that of xo an d X2 by ^a(fi'—fj,), 
but xo and X2 but remain degenerate with each other. 
We can break the remaining degeneracy using the same 
VEV of the 5-plet as in eq. (fT3|) by including the Yukawa 
interaction 

h^abXaXb (14) 

This splits X0.2 by ±/iS. The resulting spectrum has the 
form 

M 2 \ ( hE- |a(/i' - fi) \ 

Mi = M - an + (15) 

M j \ -hE - \ol{^ -11) J 

We should choose hE > m e to allow the decay of X2 — ► 
Xoe + e~, and \a(ji — fj,') = rn e — SM to get the small 
splitting SM. Using the value a = 1.7 x 10 -3 suggested 
by our analysis of the INTEGRAL signal, and assuming 
that A = E/2 for example, we get the required spectrum 
with E 11.4 GeV, A 5.7 GeV, /i~4.5x 10" 5 . 



4 It is important to put the external DM states on shell to get the 
correct result. 
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We can easily construct a potential for the Higgs sector 
which leads to the desired symmetry breaking pattern: 



V = A 1 (±£ afc £ Qb -£ 2 ) 2 + A 2 (A Q A Q - A 2 ) 2 
+ A 3 A a £ afc S bc A c (16) 



To see that this can work, first consider the limit A3 = 0. 
The VEV of A can always be rotated into the Ai direc- 
tion by a global SU(2) transformation, while the VEV 
of £ has no preferred orientation. When A 3 is turned 
on, (£ ob ) prefers to have the elements in the first row 
and column vanish, (£i a ) = (S n i) = 0, so that the A3 
term remains zero. We can still perform a global rota- 
tion around the 1 axis to make (£02) = (£20) vanish; this 
rotation leaves (A ) invariant. 

To summarize, this model is extremely simple: it needs 
a new SU(2) (not SU(2)xU(l) [29]) in the dark sec- 
tor, broken by a new Higgs triplet and 5-plet, and a 
dimension-5 coupling which induces mixing between the 
photon and the new B' gauge boson. These ingredients 
easily give us the desired mass spectrum, fig. QJa) , and 
coupling of one of the gauge bosons B' to leptons. The 
stability of the middle state \i is guaranteed by an un- 
broken Z2 symmetry, which also keeps the 0.1 GeV-scale 
B and B" gauge bosons stable, but these (as we will show 
presently) are cosmologically harmless and not subject to 
accelerator constraints. 



VI. COSMOLOGICAL IMPLICATIONS 

It is interesting to note that the model we have put 
forward based primarily on the INTEGRAL and ATIC 
observations happens to predict a relic density for the 
DM which is not so far from the required value. To match 
the WMAP value Q x h 2 = 0.1, one needs ((Tann«rel) — 1 
pb-c [30(; for three colors of nearly degenerate DM, this 
becomes 3 pb-c. Take v Ie ] — 2p\/M, where p\ is the 
momentum of one of the incoming \i ' s i n the- center of 
mass frame. At the freeze-out temperature T ~ M/20, 
pi ~ y/3TM. The Mandelstam variable t ranges between 
the values t± = ~M 2 ± 2Mp x = —M 2 ± At/2, while 
s = 4M 2 . A somewhat detailed computation gives us 
an estimate for the matrix element for the annihilation 
process shown in fig. Ufa) as \A4\ 2 = 9g 4 (however we 
have not been careful enough to determine whether the 
three diagrams interfere constructively or destructively — 
this estimate assumes the former). Then 



CTW ro l 



|A^| 2 At _ 97ra 2 



32ns Pl M 2M 2 



(17) 



For the values we have favored thus far, a = 1.7 x 10~ 3 
and M = 500 Gev, (fT?)) comes out too small by a factor 
of ~ 100. We should thus hope to find an example of 
sufficient positron production for INTEGRAL at a value 
of a approximately 10 times higher, a = 0.02. This is 
preferable to achieving the right relic density by making 





FIG. 8: Left (a): diagrams determining the relic density of 
the DM particles \i- Not shown is the w-channel version of 
the first diagram. Right (b): diagrams determining the relic 
density of the new gauge bosons Bi . 



M smaller, since that would force us to give up our ex- 
planation for ATIC. As we mentioned before, the numer- 
ical scattering computation becomes prohibitively slow 
at these larger values of a; work is in progress to explore 
this region of parameter space. 

However, it is possible to add extra interactions to the 
model to adjust the relic density without the need for 
changing our preferred value of the gauge coupling or 
DM mass. If there is an extra Z' gauge boson with the 
couplings g x Xa^7 M 75Xa to the DM and g e e R Z^e R 
to right-handed electrons, it can provide extra annihi- 
lation channels which easily bring (a ann u re i) up to the 
required value; we just need g x g e jM\, = 10g 2 /M 2 . If 
the Z' for some reason couples mainly to e R and not 
other SM particles, this can maintain the preference for 
annihilation into e + e~ but not heavier charged particles, 
as indicated by PAMELA/ ATIC. Such a Z' would have 
to correspond to a U(l) symmetry broken at the scale M 
since the Xi masses do not conserve the current to which 
Z' couples. If such a Z' solves the relic density problem, 
then its contribution to annihilation in the galaxy will 
also dominate the PAMELA/ ATIC signals, giving nearly 
monoenergetic e + e _ via the process of fig. [9j 





%0 Xl 



FIG. 9: Processes which can give a sharper feature in the 
electron spectrum, as suggested by the ATIC data. 
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The analogous freeze-out computation can be done for 
the diagrams of fig. [5](b) to determine the relic density 
of the fi = 100 MeV stable gauge bosons B-B" . The 
estimate for (cann^rd) is similar to to eq. (TT7|) , except 
that the relevant mass scale is fj, rather than M. This 
makes the cross section larger by a factor of 10 7 . Since 
f2 scales like 1/cx, the cosmological energy density of the 
stable B's is negligible, tt B ~ 10~ 6 . They can annihilate 
into e + e~ because of their coupling to B', but this would 
probably be very difficult to detect due to the low density 
of the B's. The stable gauge bosons thus seem to be 
rather innocuous cosmological relics. 

Note added. After completing the first version of this 
work, M. Pospelov pointed out a potentially serious diffi- 
culty: the process XiXi ~* XoXo can also proceed through 
Sommerfeld-enhanced scatterings in the early universe, 
exponentially depleting the density of xi states which 
we assumed to be as numerous as the xo's. This would 
take place well below the initial freeze-out temperature, 
when the kinetic energy of the %'s becomes comparable 
to the mass splitting ~ 2m e , and their velocity is of or- 
der y / 2m e /Mi ~ 10~ 3 , close to the value in the galaxy 
today. 

To estimate the severity of the problem, we compute 
the rate of X1X1 ~> XoXo scatterings over the Hubble 
rate as a function of temperature. We take the value of 
the cross section required for matching the INTEGRAL 
observation, (av) = 10~ 28 (Afi/TeV) 2 (u /v) cm 2 (where 
vq is the typical velocity in the galaxy), and the density 
n = vT 2 with v — 10 -12 x (500 GeV /Mi) corresponding 
to the value 0.dm — 0.1. The DM fell out of equilibrium 
at T = M/20 = 25 GeV; at this time its momentum 
was p ~ V3MT = My/3/20, and thereafter it redshifted 
like T. We can write p = c p T, where c p = \/60 if the 
DM underwent standard freeze-out. However we will be 
interested in more general, nonthermal ways of generat- 
ing the DM, so we keep c p unspecified for now. Taking 
H = y/g;T 2 /M p , we find 



n{av) 
H 



= 4 x 10" 



vM x M p 
g^CpGeV 2 



2 x 10 5 



(18) 



This is independent of T except through the number of 
species g*. From (fig]) it is clear that the xi -depleting 
process will be in equilibrium in the standard scenario, 
where c p = \/60. We need c p ■> 10 5 to avoid the deple- 
tion of xi- This could happen if the x's were produced 
nonthermally. A straightforward way would be through 
the late decays of a heavy scalar S — > XiXi- For exam- 
ple if trig >, 2M\ and the decay occurs at T < 5 MeV, 
then p ~ Mi and c p = Mi/T ~ 10 5 . This nonstandard 
production mechanism requires that the standard one be 
subdominant, i.e., we need larger values of a or of the Z' 
couplings (as discussed above) to suppress the relic den- 
sity at freeze-out of the perturbative annihilations. It is 
straightforward to verify that the x's will not reestablish 
kinetic equilibrium with the SM particles after this point. 
By design, the strongest interaction with the SM is the 



scattering Xi e ~* Xi e y ia ^-channel exchange of the Z' 
boson. The rate for this is 8 orders of magnitude smaller 
than the Hubble rate at T ~ 5 MeV. 

One question which arises with respect to out-of- 
cquilibrium production of DM through decays is how to 
get a sufficiently small coupling between S and x so that 
the decay happens at such low temperatures. However, 
this can be naturally explained in the context of a Grand 
Unified Theory (GUT), with a heavy gauge boson X that 
couples to S and to x- Integrating out the X boson gives 
rise to an interaction of the form g 2 x \S\ 2 x a x a /M x , where 
gx is the GUT gauge coupling. This interaction does not 
allow for S to decay, but if we suppose that S gets a VEV 
at the TeV scale, a coupling which leads to decays is gen- 
erated, A = g 2 x (S)/M X - The decay rate T = 3A 2 M s /16tt 
must equal the Hubble rate H = g r T 2 /M p at T < 5 
MeV, leading to A < 10" 12 for M s - 1 TeV. Taking 
(5) - 10 TeV, we find M x /g 2 x - 10 16 GeV, the GUT 
scale. Of course S would need to have stronger couplings 
to other particles to freeze out at the right density. 



VII. CONCLUSIONS 

We have demonstrated a modification of the original 
XDM scenario which makes the effect strong enough 
to explain the 511 keV anomaly, without any radical 
change in our understanding of the DM distribution in 
the galaxy, and within a theoretically sound class of par- 
ticle physics models. The key idea is to have two nearly 
degenerate states of stable DM, such that the heavier one 
must first be excited by the small energy SM ~ 100 keV 
into a third state, in order to decay to the lowest state 
(fig. [Tj). Moreover, we need at least two kinds of new 0.1 
GeV-scale gauge bosons, B^ and B'^ (with a theoretical 
preference for a third, B'^), such that B^ and B" mediate 
the excitation XiXi ~~ ¥ X2X2, while B'^ is responsible for 
the decay xi ~ * Xa e+e ~ ■ A novelty of this scenario is 
that the first excited DM state is stable, in addition to 
the ground state. 

Although our primary interest was to address the IN- 
TEGRAL observation, we have also been motivated by 
the suggestion [l[ of a class of models which can simul- 
taneously explain the other anomalies; thus we have fo- 
cused on parameters where M\ is at least 500 GeV. Since 
the INTEGRAL signal scales like Mj~ 4 , we can more eas- 
ily accommodate it by relaxing this requirement. To 
make contact with ATIC/PPB-BETS, it may be desir- 
able to find a model in which M\ is somewhat higher 
than the 500 GeV value we were able to obtain thus far, 
since the ATIC data indicate an excess going up to 800 
GeV. The fact that we did not yet find an example with 
Mi > 500 GeV may just be due to the numerical chal- 
lenges of our method, where the integration time to solve 
eq. ([3]) (subject to the appropriate boundary conditions) 
becomes increasingly long as T (eq. (j6])) and Mi are in- 
creased. It would be desirable to find an approximate 
analytic solution in the large T regime to make further 
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progress. 

Our proposal for modifying the excited dark matter 
spectrum has also motivated us to construct a model 
which is simpler than others that have been suggested 
along these lines [29(; for example the new gauge group 
in the dark sector is SU(2) rather than SU(2)xU(l). A 
bonus of our model is the possibility of a new signal due 
to annihilation of DM states into a single e + e~ pair in- 
stead of two pairs, giving a sharper feature in the lepton 
spectrum at energies near the DM mass. In this paper 
we have not tried to predict the PAMELA / ATIC /PPB- 
BETS signals in as much detail as we have treated that 
of INTEGRAL; doing so would clearly be a worthwhile 



next step. 
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